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ULTRASONIC ELASTOGRAPHY PROVIDING 
AXIAL, ORTHOGONAL, AND SHEAR STRAIN 



CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] This application is a continuation-in-part of U.S. application Serial No. 
10/765,293 filed January 27, 2004. 

STATEMENT REGARDING FEDERALLY 
SPONSORED RESEARCH OR DEVELOPMENT 

[0002] This invention was made with United States government support 
awarded by the following agencies: 
NIH CA86278 

The United States has certain rights in this invention. 

BACKGROUND OF THE INVENTION 

[0003] The present invention relates to ultrasonic medical imaging and in 
particular to an apparatus and method for making ultrasonic elastography 
measurements. 

[0004] Elastography is a new imaging modality that reveals the stiffness 
properties of tissue, for example, axial strain, lateral strain, Poisson's Ratio, Young's 
Modulus, and other common strain and strain-related measurements. The strain 
measurements may be collected over an area and compiled as a two-dimensional 
array of data, which may then be mapped to a gray or color scale to form a strain 
"image". Analogously, strain measurements may be collected over a volume 
displayed either three-dimensionally or as a series of stacked two-dimensional 
images. 

[0005] In quasi-static elastography, two images of the tissue may be obtained by 
the ultrasound device in two different states of compression, for example, no 
compression and a given positive compression. The tissue may be compressed by 
an external agency such as a probe or the like, or may be compressed by muscular 
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action or the movement of adjacent organs. Strain may be deduced from these two 
images by computing gradients of the relative local shifts or displacements in the 
images along the compression axis. Quasi-static elastography is analogous to a 
physician's palpation of tissue in which the physician determines stiffness by 
pressing the tissue and detecting the amount the tissue yields under this pressure. 
[0006] In dynamic elastography, a low-frequency vibration is applied to the 
tissue and the tissue vibrations are measured, for example, using Doppler detection. 
[0007] Typically quasi-static elastography produces a strain measurement only 
along the axis of compression. However, lateral strain or elevational strain (both 
perpendicular to the axial strain) may be of value both in deducing qualities like 
Poisson's Ratio and in countering the effects of lateral motion in de-correlating the 
axial displacement of the tissue. In addition, shear strain images can also be 
obtained. 

[0008] A number of different methods have been used to obtain lateral strain in 
quasi-static elastography. Some of these techniques assume knowledge about the 
compressibility of the tissue (for example, Poisson's ratio), and thus cannot be used 
for measurements in which such tissue properties are to be determined. U.S. Patent 
6,270,459 to Ophir describes a technique which interpolates between successive 
axial rays or echo signals to provide a basis for horizontal displacement 
measurement using a correlation technique. 

[0009] U.S. patent application 10/765,293 by the present inventors, entitled: 
"Ultrasonic Elastography With Angular Compounding 1 ' and hereby incorporated by 
reference, provides a method of obtaining both axial and lateral strain using multiple 
angles of ultrasonic measurement and an angle-dependent weighting factor based on 
an assumed value of the Poisson's ratio The compounding of the measurements 
from multiple angles improves the accuracy of the strain determinations. 



SUMMARY OF THE INVENTION 

[0010] The present invention provides an improved method of extracting axial, 
lateral, and elevational strain (more generally orthogonal strain tensor components) 
and displacements by developing a model representing projected displacements 
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measured by ultrasound at multiple angles and as a function of axial and lateral 
displacement. The model is then fit to actual ultrasonic measurements and a noise- 
reduced estimation of axial and orthogonal displacements are extracted from the 
adjusted model. Axial and orthogonal displacements are then used to develop other 
parameters including axial and orthogonal strain, Poisson's ratio, and shear strains. 
[0011] Importantly, the model does not need to presuppose the mechanical 
properties of the tissue (e.g., Poisson's ratio), and thus may be used in measurements 
of such mechanical properties. The fitting of the model and data provides a 
sophisticated method of reducing the effect of noise on any one measured value. 
[0012] Specifically then, the present invention provides an ultrasonic 
elastography system including an ultrasonic acquisition assembly adapted to provide 
a set of ultrasonic signals from a plurality of voxels in a region of interest at a 
plurality of angles through the voxels, the set of ultrasonic signals including a first 
subset of ultrasonic signals taken with tissue of the region of interest in a first axial 
compressive state and a corresponding second subset of ultrasonic signals taken with 
tissue of the region of interest in a second axial compressive state. A processor 
receiving the set of ultrasonic signals, executes a stored program to: (i) measure the 
displacement of each voxel projected along the angle of each of the ultrasonic 
signals between the first and second compressive states; (ii) fit a model providing 
projected displacement as a function of ultrasonic signal angle and axial and 
orthogonal displacement to the measured displacements at each voxel; and (iii) 
determine axial and orthogonal displacement for the voxels from the fit model. 
[0013] Thus one object of at least one embodiment of the invention is to provide 
an improved method and apparatus for extracting axial and orthogonal 
displacements from a multi-angled ultrasonic acquisition. The fitting of a model to 
multiple measurements reduces the effect of noise in any single measurement. 
[0014] The processor may further determine parameters for the voxels related to 
the determined axial and orthogonal displacements. 

[0015] It is yet another object of at least one embodiment of the invention to 
provide information about displacement related parameters such as Poisson's ratio, 
normal and sheer strains and stress distributions. 
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[0016] The model need not presuppose material properties of the voxels. 

[0017] It is yet another object of one embodiment of the invention to provide a 

method that can be used to determine material properties of the voxels. 

[0018] The plurality of angles of ultrasonic signals are in multiple perpendicular 

planes. 

[0019] It is yet another object of at least one embodiment of the invention to 
provide for 3-dimensional displacement measurements and parameters derived from 
such measurements. 

[0020] The ultrasonic acquisition assembly may include a transducer selected 
from the group consisting of: a single transducer element moved in location and 
angle, a multi-element transducer moved in location and angle, a phased array 
transducer sweeping in angle and moved in location, and beam steering on 
multielement transducers held stationary or moved in location. 
[0021] It is yet another object of at least one embodiment of the invention to 
provide a method that flexibly works with a variety of different transducer designs. 
[0022] The ultrasound elastography system may include a display device and the 
processor may provide an image output to the display based on the determined axial 
and orthogonal displacements, for example, an image based on axial and/or lateral 
strain, and/or elevational strain, and/or Poisson's ratio, and/or shear strain. 
[0023] It is yet another object of at least one embodiment of the invention to 
provide a system that produces new views of tissue not normally obtained from 
ultrasonic equipment. 

[0024] The foregoing objects and advantages may not apply to all embodiments 
of the inventions and are not intended to define the scope of the invention, for which 
purpose claims are provided. In the following description, reference is made to the 
accompanying drawings, which form a part hereof, and in which there is shown by 
way of illustration, a preferred embodiment of the invention. Such embodiment also 
does not define the scope of the invention and reference must be made therefore to 
the claims for this purpose. 
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BRIEF DESCRIPTION OF THE FIGURES 

[0025] Fig. 1 is a simplified block diagram of an ultrasound machine suitable for 
acquiring multi-angle ultrasonic measurements used in the present invention; 
[0026] Fig. 2 is a front elevational view of a phased array transducer that may be 
part of the ultrasound machine of Fig. 1 to provide a steerable and focusable beam to 
obtain multi-angle measurements; 

[0027] Fig. 3 is a view of the transducer of Fig. 2 as moved laterally in 
translation to provide for multi-angle measurements of a plurality of voxels; 
[0028] Fig. 4 is a simplified representation of three-dimensional ultrasonic 
measurements obtained by moving the transducer Fig. 2 both laterally and in 
elevation in two perpendicular sweeps to obtain multiple angled measurements in 
two perpendicular planes; 

[0029] Fig. 5 is a geometric diagram of displacement of tissue in a voxel 
showing various quantities described in the present invention; 

[0030] Fig. 6 is a plot showing two curves of a model of projected displacement 
for different displacements as a function of angle superimposed with measured data 
points; 

[0031] Fig. 7 is a flowchart showing steps of a program executed by a processor 
of the ultrasound machine of Fig. 1 in providing axial and orthogonal displacements 
for each voxel and parameters based on those displacements; and 
[0032] Fig. 8 is a geometric diagram showing quantities used in the extension of 
the techniques of Figs. 5 and 6 to three dimensions. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

[0033] Referring now to Fig. 1, an elastographic imaging system 10 suitable for 
use with the present invention may employ acquisition electronics 12, for example, a 
commercial ultrasonic imaging machine programmed according to the present 
invention, or a combination of such a machine and an external computer processing 
data acquired by the commercial machine. 

[0034] Generally the acquisition electronics 12 includes an ultrasonic transducer 
18, interface circuitry 20, and a processor 16 communicating with a terminal or 
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display device 22. In operation, the transducer 18 provides ultrasound signals and 
position signals to the interface circuitry 20 which provides digitized versions of 
both suitable for processing by processor 16 to display an image on display device 
22. 

[0035] Together, the acquisition electronics 12 provides the necessary hardware 
and software to collect a series of ultrasonic echo signals at a plurality of angles 14 
through an imaged object 24, such as an organ. Two acquisitions are acquired at 
each angle 14 with different degrees of compression of the imaged object 24 in an 
axial direction by a compressor 26. Compression amounts are typically about 1% 
for tissue and may, for example, be an uncompressed and compressed state or two 
compressed states with different amounts of compression. 
[0036] As used herein "axial" will be the direction that is parallel to the 
compressive force provided by compressor 26 along a z-axis and the term 
"orthogonal" will be used to describe either or both of a "lateral" direction 
perpendicular to the axial direction along an x-axis and "elevational" direction 
perpendicular to the axial and lateral direction along a y-axis. As will be understood 
from the following description the present invention is applicable to both lateral and 
elevational strains. 

[0037] Under compression, the imaged object changes by decreasing in axial 
height by a decrease amount 27 and expanding laterally by expansion amount 28 
determined generally by a Poisson's ratio for the material. The Poisson f s ratio may 
allow for differentiation between normal and abnormal tissues. Further, for 
poroelastic materials, the Poisson's ratio may change during compression because of 
unbound liquids leaving the material, and thus imaging of Poisson ! s ratio may be 
used for quantitative assessment of fluid transport in regions of edema, 
inflammation, or the like. 

[0038] Referring now to Fig. 2 in a preferred embodiment, the transducer 18 is a 
multielement array transducer having multiple transmitting and receiving elements 
32 such as may create independent wavefronts 34 as will be understood to those of 
ordinary skill in the art. A group of these elements is activated to create an 
individual ultrasound beam, and the beam's central axis is translated along the array 
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to create an ultrasound image. Control of the phase of wavefronts of the elements 32 
can provide for a steerable beam 33 that allows multiple angles of acoustic signals to 
be transmitted from and received by the transducer 1 8 without movement of the 
transducer 18. The central axes of these steered beams may also be translated along 
the array to create an angled image. Thus, signals from each voxel are obtained from 
a plurality of beam angles. 

[0039] Referring now to Fig. 3, moving the transducer 18 lateral direction (to a 
position denoted by transducer 18') allows individual voxels 30 to be measured at 
multiple angles 14 determined by the number of lateral movements and acquisitions. 
For example, each voxel may be measured by beams ranging from +45 to -45 
degrees in one-degree or other fixed angular increments. 
[0040] Referring now to Fig. 4, a set of scans with incremental lateral 
displacements along the x-axis may produce a scan plane 35 suitable for a two- 
dimensional analysis as will be described below. Moving the transducer 1 8 without 
rotation about the axial direction, after each scan plane 35 is acquired, allows the 
acquisition of multiple scan planes 25 covering a volume of the imaged object 24. 
In each of these scan planes 35, axial and lateral displacements can be determined. 
Rotation of the transducer 18 by 90° with incremental elevational displacements 
along the y-axis may produce a scan plane 37 perpendicular to scan planes 35. 
Moving the transducer 1 8 without rotation about the axial direction, after each scan 
plane 37 is acquired, allows the acquisition of multiple scan planes 37 covering the 
same volume of the imaged object 24 as the scan planes 35, but allowing axial and 
elevational displacements to be determined. 

[0041] Instrumentation attached to the transducer 18 or analyses of the 
ultrasonic signals acquired during the scans may be used to provide an indication of 
the angle 14 and relative location of each acquired ultrasonic signal of the scan 
planes 35 and 37. Additional details with respect to an elastographic imaging 
machine suitable for programming per the present invention is described in co- 
pending U.S. application cited above. 

[0042] Referring now to Fig. 5, for each given voxel 30 within the imaged object 
24 the axial compression will cause a displacement both in axial and orthogonal 
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directions so that the volume element 30 under compression moves to a position 
denoted by volume element 30 f typically having an axial or z-component portion d z 
and a lateral or elevational component d x or d . When this displacement denoted 
by vector d is measured along an angle 14 of one ultrasonic signal, only a projection 
of the displacement d will be detected by that angle 14. For small voxels 30, d 
may be approximated by a straight line and this projected displacement q e will be as 
follows: 

q e = d 2 cos 0 + d x sin 0 + n Q (1) 

where d z and d x are axial and lateral components of the displacement vector 
d in the z- and x- directions, respectively, 0 is the angle 14 with respect to the z 
axis, and n 0 is noise or inaccuracy in the measurement process. Minor 
modifications of equation (1) will provide the relationship for axial and elevational 
displacement. Note that this equation (1) does not require knowledge of Poisson's 
ratio or any other mechanical property of the imaged object 24 but is simply the 
result of geometric principles and is thus material independent. 
[0043] Referring now to Fig. 6, a family of curves 40, 40* (only two are shown 
for clarity) may describe noiseless measurement of an arbitrary displacement vector 
d as a function of 0. Each curve has a generally downward opening arc and the 
following the equation: 

Pe - d 2 cos 0 + d x sin 0 (2) 
[0044] Generally as will not be described, in the present invention, for each 
voxel 30, measured data along the angles 14 is used to select one of these curves that 
best fits the measured data and this curve and its equation reveal an estimate of 
noiseless measurement of d z and d x (or d ) 

[0045] Referring to Fig. 7, in a first step of the present invention indicated by 
process block 42 such as may be implemented by software running on processor 16, 
a pair of ultrasonic signals s 0 and s^are acquired along each of different unit 
vectors u e at angles 0 during a multi-angle echo signal acquisition. Signal s e is 
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obtained under a first state of tissue compression (e.g., no compression) and signal 
s' 0 is obtained with a predetermined compression. 
[0046] These signals s 0 and s f d are processed according to standard 
elastographic techniques (for example, cross-correlation) to deduce a projected 
displacement g^at each voxel 30 for each angle 0 per equation (1) above. 
Generally, the ultrasonic measurements along different angles 14 will not 
necessarily cross the center of the voxels 30 and accordingly interpolation may be 
applied to the displacements obtained at different angles 14 to bring them to a 
common point in the voxels 30. 

[0047] Per process block 46, for each voxel 30, the measurements q 0 are fit to 
one of the family of curves providing values of p e that are closest to the values of 
q Q . While this may be done, for example, by iterating through combinations of d z 
and d x for the range of angles 0 , however, in the preferred embodiment, a least 
squares fit is rapidly obtained by using the following matrix relationship: 

d=(A T A)~ X A T q (3) 
where: 



d is the displacement vector 



q is the set of measured projections of displacement 



4* 



and 



cosO x 


sin 9 X 


cos0 2 


sin 0 2 


cosO m 


sin 0 m 



A = 



[0048] The matrix product A 7 A is called the Grammian or Gram matrix and a 
solution exists for this equation (3) because the columns of A are linearly 
independent. 
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[0049] At process block 47 the axial and lateral displacements d z and d x are 
taken from the solution of equation (3) and then, as indicated at process block 48 
used to deduce, for example, axial and lateral sheer ( e 22 , e xr ) being change in 
displacement as a function of the z or x axis, respectively, per the following 
equations (4) and (5): 

dx 

[0050] In this case an assumption is made that the stress field is uniform, 
however, these equations may be easily modified as will be understood to one of 
ordinary skill in the art to accommodate varying stress fields. 
[0051] From these quantities Poisson's ratio, shear strain (e zx ), and other 
parameters may be readily deduced as will be understood to those of ordinary skill 
in the art. Shear strain is given by the following equation: 

lfdd 2 dd x } 

[0052] Poisson's ratio may be determined for each voxel by division of the strain 
in the x direction and by the strains in the z direction recognizing that this method 
provides a good approximation only if the voxels are subject locally only to axial 
compression. Generally this will not be true around inclusions and therefore the 
present technique can be used to emphasize such features. 

[0053] Each of these deduced parameters can be used to produce an image, for 
example, a Poisson's ratio image on display device 22. 

[0054] It should be noted that the present invention does not preclude the use of 
other models than that provided by equations (1) and (2), including, for example, 
those that take into account other features of the measurement process. Note that the 
model may accommodate different numbers of ultrasonic measurements and 
accordingly will work even when the angles 14 along which measurements may be 
obtained are limited by structure. 
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[0055] Referring now to Fig. 8, the two dimensional models described above 
may be expanded to three dimensions by analyzing projected displacement p 9<f> as a 

function of both 0 in a lateral plane and <f> in an elevational plane perpendicular to 
the lateral plane as shown in Fig. 8. In this case curves 40 through 40' representing 
the modeling will be surfaces per the following equation: 

Pe,<f> — d 2 cos 0 cos <fi + d x sin 0 + d cos 0 sin (j> . (7) 
[0056] The fitting process takes measurements q e4 and fits them to one of these 
curves using equation (3) with the elements defined as follows: 



A = 



d = 



4*4 



cos 0 } sin (j> x sin 0 X cos 0 X sin (f> x 

cos 0 X sin <f> n sin 0, cos 6 X sin <j> n 

cos 0 2 sin <f> x sin 0 2 cos 0 2 sin (f> x 

cos e m sin <f> n sin 0 m cos 0 m sin ^„ 



[0057] It will be apparent from this description to one of ordinary skill in the art 
that a number of variations may be made in the preferred embodiment that will still 
be within the spirit and scope of the claims. For example, while the preferred 
embodiment contemplates ultrasonic echo signals, the present invention can also be 
used with through transmitted ultrasonic waves. The invention contemplates use 
with alternative ultrasonic transducer technology and other methods of obtaining 
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multiple angles of ultrasonic measurements through the voxels of the imaged object 
24 including area array transducers providing cone beams that may be phased to 
sweep through angles over three dimensions or individual ultrasonic transducers that 
may be moved in translation and rotation to acquire the necessary data. 
[0058] It is specifically intended that the present invention not be limited to the 
embodiments and illustrations contained herein, but include modified forms of those 
embodiments including portions of the embodiments and combinations of elements 
of different embodiments as come within the scope of the following claims. 
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